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SUMMARY 


The  electromagnetic  field  strength  measure  nents  ) 

obtained  during  the  1966  FCC/USAl’  POPSI  Project  have  been  classified 
according  to  propagation  mode  and  the  data  from  periods  of  precipita¬ 
tion  have  been  analyzed  in  detail.  Cumulative  probability  distribu¬ 
tions  were  generated  for  the  effective  radar  reflectivities  derived 
from  the  bistatic  electromagnetic  measurements  and  compared  with  the 
probability  distributions  of  the  surface  rainfall  rates  derived  from 
the  accumulations  of  the  United  States  Weather  Bureau  recording  rain 
gauges  in  the  area.  The  distribution  functions  were  then  adjusted 
by  means  of  a  least  squares  regression  line^^to  obtain  a  Z(p) 
relationship  for  effective  reflectivities  i^  excess  of  10°  mm  °/rir 
and  surface  rainfall  rates  up  to  105  mm/hr.  ^.The  relationship  thus 
obtained  has  been  compared  with  other  Z-R  relationships  based  upon 
the  analysis  of  drop  size  distributions  and  has  been  tested  against 
Independent  rain  gauge  data  in  the  POPSI  Project  area.  The  final 
approximation  resulted  in  a  standard  deviation  for  estimating  Z(p) 
from  R(p)  of  less  than  1.6  dB  for  the  New  Jersey  rainfall  data^Nand 
appears  to  be  more  representative  than  Z  for  the-actijal 

relationship  between  Z(p)  and  R(p)  for  convective  storms,  ^he  alti¬ 
tude  dependence  of  the  reflectivity  from  precipitation-connected 
phenomena  in  the  New  Jersey  coastal  area  has  been  demonstrated  and 
discussed  to  some  extent. 
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1.0  INTROPUCTION 

During  the  period  froo  15  February  1966  to  16  Februarx*^  1967. 
the  FCC  aud  USAF,  with  the  coopa-^ation  or  the  0.  3.  Coast  Guard, 
the  FAA,  the  U.  S.  Weather  Bureau,  and  NASA,  conducted  an  invest! 
gatlon  of  the  signal  power  scattered  fro-ji  precipitation  and  other 
mechanisms  iu  the  common  voluines  established  by  the  inrerssctions 
of  the  beams  from  transmitting  auceana  op^reting  in  a  configura- 
tiou  sinulatlng  a  satellite  oarth  stction  aod  receiving  catennas 
configured  in  i  manner  cyp:.tal.  of  tc-rtestrial  rBlcro::ave  radici- relay 
stations.  The  investigation  was  conducted  at  a  rr-diation  wavelength 
of  5.?1  centimeters  (5.75  GtU)  in  &u  area  near  the  New  Jersey  Coast 
and  was  designed  to  obtain  da;::a  ror  a  statistical  tteatmeat  of  the 
scattered  interference  prcblein.  The  details  of  this  ^recipi* 

tation  and  Off  Path  ^c«ctered  ^nserferent^ )  Project  may  be  found  i? 
FCC  Research  Division  Report  Ko.  R-6S0i,  dated  15  Karclv  1968.(1) 

lu  the  period  immediately  following  the  I'CC/USAF  project, 
several  domestic  and  international  groups  were  conyenca  for  the 
purpose  of  agreeing  upon  the  form  and  direction  of  research  pro¬ 
jects  to  further  the  investigation  of  the  off-^path  ptopagaticn 
phearnaena  acd  to  derive  allocfvtion  criteria  which  would  recog¬ 
nize  the  interference  potential  to  terrestrial  microwave  aud 
satellite  consnuai cation  system  earth  stations.  However,  it  was 
not  untxl  midway  through  fiscal  year  1970  that  a  project  was 
Initiated  uc,de:r  the  management  of  a  U.  S.  interagency  group  (2), 
Meanwhile,  the  FCC  Research  Division  decided  to  proceed  with  a 
detailed  analysis  of  the  original  POPSI  Project  data, 

2 . 0  SEPARATION  AND  GROUPING  OF  DATA 

Although  the  original  objectives  of  the  FOKI  Project  did 
not  include  the  positive  identification  of  the  dominant  propagation 
mechanism  at  all  times,  the  raw  data  records  and  meteorological 
data  inputs  were  adequate  to  permit  a  computerized  format  for  the 
identification  of  the  predominant  propagation  mode  during  certain 
specific  time  periods. 

Without  going  into  detail,  the  tools  used  for  propagation 
mode  identification  included: 

a.  Comparison  of  the  medians  of  simultaneous  5-minute  Intervals 
of  the  great-circle  and  off-path  signals,  and  correlation  of 
the  differences  with  known  antenna  side- lobe  radiation  patterns. 

b.  Chart  recordings  from  the  U.  S.  Weather  Bureau  rain  gauges  in 
the  vicinity  of  the  propagation  paths. 
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c.  V^sttccil  profiles  of  hu333idity»  temperature  aAfI  refractive 
iadex  gradieacs  constructed  from  radiosonde  data  taken  every 
four  hours  from  a  location,  neat  the  transmitter'  site  aad  Intel 
mittfiutly  rroa  a  location  oa  the  grsaC-circle  path. 

d.  Time-lapse  photographs  of  the  ?PI  scope  of  a  57  weather 

radar  located  on  the  great-circle  path.  (5>  >»  5) 

e.  Correlation  of  the  short-term  variability  {fading  rate)  of 
the  cff-patb  signal  with  that  of  the  great^circle  path 

signal,  no,  15,  5A) 

For  the  detailed  analyses,  the  data  were  divided  into  hourly 
segments  and  assigned  to  one  of  three  groups,  according  to  the 
propagation  mode. 

1.  Precipitation  Scattering 

1.  Guided  Propagation  (ground-based  and  elevated  layer 
ducting) . 

3.  Mixed  (coebination  of  precipitation  scattering  and  guided 
propagetioii) . 


Very  few  difficulties  were  encountered  in  the  Identification 
and  analysis  of  off-path  signal  enhancemonis  attributed  to  precipi¬ 
tation  alone,  and  the  periods  of  ground-based  ducting  were  easily 
recognized.  However,  set^ere  problems  In  Ideutiflcatlon  and  analyses 
i^ere  experienced  with  the  high  signal  levels  associated  with  elevated 
temperature  inversions  or  humidity  lapses  in  the  abcccce  of  measur¬ 
able  surface  precipitation  and  with  the  direcrional,  part tally- coher¬ 
ent,  propagation  often  associated  with  precipitation  from  s  heavily 
stratified  troposphere. The  analysis  of  this  'mixed  mode" 
propagation  has  been  difficult  because  of  insufficient  metaorologicel 
input,  unfavorable  path  geometry,  and  the  lack  of  an  adequote  model. 


3.0  COSTEMPORARY  rilEORY 

The  scattering  and  attenuation  of  electromagnetic  wavia  by 
particles  in  the  atmosphere  are  complex  functions  of  the  particle 
size,  dielectric  properties,  and  the  radiation  wavelength. 

The  theoretical  treatment  of  the  relationships  among  forworc- scattered 
power,  back-scattered  power,  radar  reflectivity,  back-scar teriag 
cross-section,  and  rainfall  rate,  are  usually  simplified  by  beginning 
with  the  Rayleigh  approxirsation  for  the  back- scattering  cross-section 
of  d  single  spherical  particle  (raindrop,  ice  particle,  hailstone, 
etc.)  having  a  diameter,  D.  This  cross-section  is 
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(3-1) 


whs?e  %.  i(f  »  function  of  the  particle  refractive  index,  and  X  is 
the  radiation  wavelength.  Equation  (3-1)  is  generally  considered 
to  be  valid  when  sP/X  0.2  ic  thiw  f*pecial  case,  the  back- 
scatter  is  prop /rtlooal  to  |K|  where 


2  1 

nv  -  1 

2  , 

m  -V  ^ 


(3-2) 


fi  is  the  square  root  of  the  complex  dielectric  constant,  e  ,  and 
is  eqT-al  to  n-j*  where  n  is  the  phase  refractive  index  and  «  is  the 
absorption  coefficient  o£  the  particle  substance.  Although  the 
experimentally-derived  values  of  n  and  k  have  exhibited  some  d^'pen- 
dence  upon  temperature  and  war-eleagth,  the  variations  in  the  centi- 
icetric  land  appear  to  be  slij^ht  and  |K|  ‘  is  usually  assumed  to  be 
0.93  for  water  and  0.176  for  ice. 

For  w  particles  per  xoiit  volume,  the  back-scattering  cross- 
section  per  unit  volume  (reflection  coefficient)  is 


U  «=  2<r  “ 


(3-3) 


If  all  the  particles  are  assumed  to  he  the  same  size,  the  quantity 
known  as  the  equivalent  reflectivity,  Z,  is  expressed  by 


Z  « 


(3-4) 


If  the  particles  are  not  the  same  size  but  instead  have  a  "drop 
size  distribution," 


Z  =  2  N^Di 
i=l 


(3-5) 


In  any  event,  eqtiation  (3-3)  becomes 


n  -  iT  !Ki^Z 


(3-6) 


It  can  be  seen  from  equttiens  (3-3),  (3-4),  and  (3-5)  that  the 
particle  diameter  is,  by  several  orders  of  magnitude,  the  most 
significant  parameter  in  the  determlnvition  of  Z,  and  hence  the 
determination  of  H  by  the  indirect  method  (conversion  of  surface 
rainfall  rate  to  reflectivity  factor,  Z,  or  to  reflection  cceffi- 
cienc,  ii  ) . 


Nearly  all  recent  investigators  are  agreed  that  there  appears 
to  be  no  unique  drop-size  distribution  for  a  given  rainfall  rate 
R;  therefore,  there  can  be  no  unique  relationship  between  Z  and  R, 

(17,  18,  19,  20,  21,  22,  23,  24)  However,  in  considering  the  effect 
of  precipitation-scattered  interference  upon  satellite  and  terres¬ 
trial  microwave  cooHnunications  systems,  the  uoe  of  empirical  rela¬ 
tionships  between  effective  radar  reflectivities  and  nearby  surface 
rainfall  rates  is  especially  attract:$ve,  inasmuch  as  the  probability 
distribution  functions  of  surface  rainfall  rates  are  parameters  that 
can  be  derived  from  meteorological  records  available  in  a  large  part 
of  Che  world. Most  of  the  relationships  between  effective  reflec¬ 
tivity  and  surface  rainfall  are  of  the  form 

Z  =  3R^ 

where  Z  is  the  effective  reflectivity  and  R  is  the  surface  rainfall 
rate.  The  empirical  constants  nre  a  and  b.  One  of  the  most  popular 
(sf  the  relationships  was  derived  from  the  application  of  the  Rayleigh 
particle  scattering  theory  to  the  observed  drop  size  distributions 
for  -/arious  rainfall  rates 


Z  =  200  R 
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where  Z  is  in  mm  m  and  R  is  in  mm/hr.  The  reflection  coefficient, 
,  may  also  be  determiufed,  more  or  less,  from  actual  electromagnetic 
waT'i.  propagation  measurements  by  tising  the  relationship. 


(3-8) 


where 


Py  =  power  at  receiver  input  terminals. 

Pt  =  power  at  transmitter  output  terminals. 

-  transmitting  and  receiving  antenna  gains 
(over  isotropic). 

\  -  radiation  wavelength. 

dt  =  distance  from  transmitter  site  to  scatterer. 
=  aistar.ee  from  receiver  site  to  scatterer. 


V 


common  volume 


B 


-  5  - 


Actually,  equation  (3-8)  is  a  re-arrangement  of  the  bl-st&tlc  radar 
equation,  which  is  in  turn  a  slight  modificatioa  (2  distances  In¬ 
stead  of  1;  2  antenna  gains  Instead  of  1;  and  pulse  length  Instead 
of  1/2  pulse  length)  of  the  radar  back- scatter  equation.  Ifflplicit 
in  the  use  of  (3-8)  for  a  forward-scatter  solution  is  the  assunu  tiori 
that  the  scattering  is  isotropic  or  exhibits  insignificant  anisotropy. 
It  is  a  scalar  formula  and  can  be  used  only  to  obtain  va7uec  of  b  or 
Z  if  they  are  also  scalar  functions.*  (26,  27,  28) 


POPSI  PROJECT  DETAILED  ANALYSIS  TECHNIQUE 


After  grouping  the  data  as  described  in  section  2.0  and 
eliminating  the  direct-path  signal  components,  the  medians  of  the 
five  minute  power  Inputs  recorded  by  the  off-path  receivers  were  cou¬ 
pled  with  the  system  parameters  to  calculate  the  effective  reflec¬ 
tivities  for  each  common  volume  configuration  by  using  the  equations, 
(3-8)  and  (3-6). 


1}  =  (4ff)^Pr^^t^r)^ 

P^GtCr^^V 


(3-8) 


Z  (eff.)  = 


(3-6) 


*  /  For  the  prediction  of  interference  to  services  having  high 
reliability  criteria,  the  infrequent  but  extremely  high 
signal  levels  exhibiting  directional,  partially  coherent 
characteristics  should  be  considered,  although  the  number 
of  patameters  in  the  explanatory  theories  is  almost  prohi¬ 
bitive.  Nearly  all  of  these  theories  involve  vector  approaches, 
such  as  geometrical  optics  approximations  or  quasi-exact 
solutions  of  the  wave  equations.  In  fact,  even  the  scatter¬ 
ing  theory,  when  larger  particles  are  involved  >  0.2), 
forces  exact  solutions  of  the  Mie  equations.  A  l^rge  segment 
of  the  POPSI  Project  data  fall  in  this  category  hut  the 
analysis  was  considered  to  be  beyond  the  scope  of  this  report. 


P'-  * 
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where  the  lA&rasteters  are  a*  previously  indicated  except  that,  in 
this  case, 

dj.  =•  distance  from  transmitter  oilte  to  center  of 
comon  volisae 

df  »  distance  frost  receiver  site  to  center  of  comson 
volume 

V  =>  cosiaon  volume  defined  by  antenna  beam  half-power 
points  and  pulse  length 

{X|^  -  0,93 

Equation  (3-7)  was  then  used  to  normalize  the  effective  reflectivity 
to  the  apparent  rainfall  rate  in  mm/hr.** 

Z  =  200  (3-7) 

The  rainfall  rates  were  derived  from  the  accumulations 
recorded  by  U.  S.  W,  B.  rain  gauges  at  NAFEC,  Cape  May,  Glassboro, 
Freehold,  Kightstown,  l^berton,  and  Marlton,  New  Jersey.  The  loca¬ 
tions  of  the  POPSI  Project  conmsoc  volumes  with  respect  to  the  rain 
gauges  are  shown  on  Figure  1.  The  MESCWET  stations  in  the  NAFEC 
area  are  also  identified,  although  technical  problems  dictated 
against  the  use  of  precipitation  data  from  them  in  the  final  analy¬ 
sts.  The  rain  gauge  recorder  charts  were  rather  difficult  to 
analyze  but  could  generally  be  interpreted  to  within  .01  of  an  inch 
for  each  15-minute  interval.  For  purposes  of  analysis. 

Precipitation  Rate:  R(iDm/hr.)  =  100  A  (^”1) 

where  A  =  hundredths  of  an  inch  accumulated  during  each  15-minute 

period  beginning  on  the  hour,  quarter-hour,  half-hour,  etc. 

This  rate  was  then  assumed  to  be  constant  during  the  entire  15-mlnute 
period. 

**  /Although  this  was  done  only  for  comparison  purposes,  it  is  a 
common  procedure  in  calculating  rainfall  rates  from  radar  re¬ 
turns.  However,  it  may  not  always  be  a  valid  pro¬ 

cedure  due  to  the  statistical  manner  in  which  the  Z-R  rela¬ 
tionship  was  obtained.  The  form  of  Z  =  200  R^*  implies  that 
it  was  the  best  fit  for  a  regression  line  of  the  dependent 
variable,  log  Z  upon  the  Independent  variable,  log  R.  Although 
the  variables  can  be  interchanged,  the  best  fit  regression 
line  may  not  be  the  same.  Futhermore,  the  greatest  deviation 
could  be  expected  at  the  high  reflectivity  extremes  of  the 
regression  lines.  (29) 
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5.0  CORREUTIffll  OF  Z  SITH  R  BY  l^TA  TIhE  PE5I0!» 

Point-by-point  correlacion  ,pf  ineasured  rainfall  races  vitb 
those  calculated  from  the  scatterfed  electioaagnetic  field  in£ro- 
duces  certain  requircioents.  (5^,  33,  34,  35)  In  addition  to  exact 
time-frame  synchronizatl  ;i,  the  rain  ^uge  systen  aperture  should 
l>e  equated  with  that  of  the  corason  volume  to  avoid  undesirable 
averaging.  Futhermore,  the  rain  gauge  field  should  be  properly 
located  with  respect  to  the  common  volume.  None  of  these  require¬ 
ments  are  likely  to  be  met  by  any  significant  research  program 
or  as  boimdary  conditions  for  an  allocation  model.  Figure  1  In¬ 
dicates  that  the  ?0PSI  Project  was  certainly  no  exception  in  there 
respects  although  the  recording  rain  gauge  facilities  were  as  exten¬ 
sive  as  can  reasonably  be  expected  outside  of  the  concentraced, 
special-purpose  networks  being  operfetec  in  a  few  limited  areas.  (57,36) 
Real  time  ccrrelatlon  cannot  be  expected  between  the  ground  based 
rain  gauges  and  effects  in  the  coEson  voltimes,  since  the  gauges  are 
displaced  from  the  common  voltimes  both  horizontally  and  vertically. 
Although  in  most  areas  it  can  be  expected  that  the  statistics  of 
rainfall  on  the  ground  may  be  very  sinsi liar  to  those  In  the  common 
volume,  the  rainfall  occurlng  here,  et  varying  altitudes  during  a 
scattering  event  may  bear  no  relation  to  the  rainfall  measured  at 
a  gauge  on  the  earth's  surface.  Further,  it  should  be  noted  that 
there  are  areas  of  the  country  where  *'dry"  thunderstorms  are  common. 

In  this  event  rain  or  hail  can  be  present  in  the  common  volume  with 
a  complete  absence  of  precipitation  on  the  ground. 

Although  the  possiblity  of  acceptable  short-term  Z-R  corre¬ 
lation  seemed  to  be  remote  because  of  the  ambiguities  Introduced 
by  space  and  time  differences  between  the  meteorological  and  the 
electromagnetic  data  inputs,  some  effort  was  made  in  this  direction. 

For  both  15-minute  and  one  hour  periods  the  effective  equivalent 
reflectivities  were  calculated  from  the  highest  median  power  of 
the  five  minute  samples  during  the  interval.  If  the  rainfall  rates 
derived  from  the  maximum  accumulations  in  the  rain  gauge,  system 
during  the  same  periods  were  at  least  one  mm/hr.*,  they  were  paired 
with  the  corresponding  reflectivities  to  obtain  the  scatter  diagrams 
nf  Figures  2  and  3. 

The  effective  equivalent  reflectivities  as  calculated  from 
equations  (3-6)  and  (3-8) ,  and  the  rainfall  rates  derived  from 
the  rain  gauge  accumulations,  were  assimied  to  be  related  by  the 
form,  (16) 

Z  =  aR^ 


*  /  The  minimum  requirement  of  1.0  mm/hr.  in  the  rain  gauge  system, 
reg'ardless  of  the  corresponding  reflectivities,  resulted  in  the 
removal  from  consideration  of  more  than  half  of  the  data  including 
13  of  the  highest  effective  reflectivities  (Z  >  10^). 
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vheze  a  anJ  b  are  constants.  By  the  logarithmic  transfomation 
the  form  becoces  linear  with  log  Z  as  the  dependent  variable,  and 
leg  R  as  the  independent  variable.  Since  a  ntmber  of  proposals 
have  been  advanced  fci'  predicting  "2"  values  from  "R"  values  mea¬ 
sured  at  the  surft^ce,  the  least  squares  regression  lines  of  log  Z 
upon  log  R  for  the  scatter  diagrams  are  of  interest.  These  have 
beer  calculated  and  are  shown  on  the  diagrams  along  with  the  loci 
of  other  empirical  relationships  of  the  same  form.  If  the  relation¬ 
ship  between  "Z"  and  "R"  were  nenrly  linear  or  if  the  data  points 
were  evenly  distributed  across  the  interval  of  consideration,  the 
least  squares  regression  lines  of  Figures  2  and  3  might  be  useful 
for  rough  prediction  on  a  point-by-point  basis.  However,  the  re¬ 
lationship,  Instead  of  being  linear,  is  in  the  form  of  a  geometric 
curve  with  the  bulk  of  the  points  at  the  lower  reflectivity  values 
which  are  of  little  Interest  and  where  the  bivariate  distrlbTition 
is  truncated  with  respect  to  rajLufall  rates.  This  results  in  a 
least  squares  rsgtessioc  line  which  is  heavily  biased  *'y  data  cf 
no  real  importance.  It  is  Interesting  to  note  that  the  hourly 
data  shows  significantly  better  correlation  than  the  15-minut8 
data.  This  was  expected  because  of  the  averaging  Involved.  Data 
for  periods  longer  than  one  hou-;-  would  probably  be  better  correlated 
and  the  correlation  of  instantaneous  Z/R  data  would  probably  be 
almost  nil. 


Hourly  Data 

Least  squares  fit: 
Coefficient  of  correlation: 
Standard  error  of  estimate: 
(Maximum  Z  from  maximum  R): 
Standard  error  of  estimate: 
(Using  Z  =  200  R^*^) 
Standard  error  of  estimate: 
(Using  Z  =  127.7  r2.26) 

15-minute  Data 


Least  squares  fit: 
Coefficient  of  correlation: 
Standard  error  of  estimate: 

Standard  error  of  estimate: 
(Using  Z  =  200  R^*^) 
Standard  error  of  estimate: 
(Using  Z  =  127.7  r2*26) 


Z  »  431,9  R*^^^ 
.34086 

6.5  dB 

8.2  dB 

10.5  dB 


Z  =  431.2 
.17714 
7.2  d3 


5.8  dB 
12.3  dB 


. .  ’ 
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The  futility  of  point-by-poiat  correlation  of  reflectivities  with 
nearby  surface  rainfall  rates  is  emphasized  by  Figure  4,  a  scatter 
diagram  of  the  data  for  October.  The  refractive  index  profiles 
for  this  period  indicated  very  little  low-level  atmospheric  strati¬ 
fication  and  the  surface  rainfall  rates  appeared  to  be  uniform 
over  a  large  area.  Theoretically,  this  period  should  promise  the 
best  2-R  correlation  for  heavy  rainfall  rates.  Actually  however, 
the  cc'-relatlon  was  quite  poor  with  reflectivities  generally  lower 
than  would  be  expected  from  the  corresponding  surface  rainfall  rates. 
This  may  have  been  caused  by  common  volumes  which  were  above  the 
0°C.  altitude  level  and  confined  to  an  area  that  was  quite  small 
compared  to  the  rain  gauge  system  aperture. 

6.0  2-R  RELATIONSHIP  BY  CORRELATION  OF  PROBABILITY  DISTRIBUTIONS 

Since  it  has  been  implicitly  issumed  that  Z  and  R  are  related, 
their  probability  density  functions  must  also  be  related  if  they 
have  been  determined  from  data  extending  over  the  same  time  frames 
and  confined  to  the  same  synoptic  weather  prediction  area.  This 
appears  to  be  the  best  approach  to  the  problem  of  Z-R  correla¬ 
tion,  since  the  effects  of  the  real  time  and  space  differences  be¬ 
tween  the  meteorological  and  the  electromagnetic  data  inputs  are 
minimized. 

Precipitation  rates  derived  from  the  15-minute  acciimulatlons 
of  I'ain  gauges  in  the  U«  S.  W.  6.  system  were  used  to  generate  a 
cumulative  probability  distribution.  The  data  from  all  of  the 
gauges  were  grouped  together  for  periods  during  which  the  transmitter 
and  at  least  one  off "path  receiver  were  in  operation.  This  distribu¬ 
tion  involved  over  87,000  15-iSiinute  samples  extending  from  February 
through  October,  and  includes  derived  surface  rainfall  rates  up  to 
105  mm/hr.  This  probability  distribution  is  specifically  indicated 
on  Figure  5  and  is  represented  by  a  smooth  curve  on  Figure  6.  The 
effective  reflectivities  were  calculated  in  accordance  with  section 
4.0  and  used  to  generate  probability  distributions  by  common  volume 
groups  as  indicated  in  Table  1.  These  probability  distributions 
are  shown  in  Figures  5  and  6,  In  these  Figures  the  ordinate  values 
for  Z  have  been  normalijsed  to  the  rainfall  rate,  R,  by  using  the 
relationship,  Z  *  200  R^'°,  In  the  interval,  40  <  R  <  100  (rom/hr.), 
the  distribution  function  for  the  15-minute  rainfall  rates  compares 
quite  favorably  with  what  a  recent  investigator  derived  for  the 
^.nstantaneous  rates  in  dense  rain  gauge  systems  in  the  same  area  (36). 

Upon  examination  of  Figures  5  end  6  it  is  apparent  that  the 
relationship  between  the  distribution  function  of  the  effective 
reflectivities  and  that  of  the  surface  rainfall  rates  is  somewhat 
dependent  upon  the  common  volume  altitude.  This  dependence  is 
not  significant  at  rainfall  rates  of  less  than  10  mm/hr.  or  with 
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connson  voltmes  beXcv  an  altitude  of  two  kllon&ters>  These  same 
boundary  conditions  might  be  ctusldered  to  be  descriptive  of  tha 
’-lidespread  uniform  precipitation  associated  with  frontal  or  oro¬ 
graphic  lifting  and  conanon  tc  temperate  climates  ia  early  spring 
and  late  fall.  For  this  type  of  precipitation,' the  relationship 
between  Z  and  R  or,  more  specifically,  the  zda tionshlp  between 
Z(p)  and  R(p)  can  be  appronisMted  by  the  ezpreosion;  Z  »  200 
There  ie  also  evldencve  that,  for  a  climate  similar  to  that  of  the 
New  Jersey  coast,  the  elimination  of  cocanon  volumes  in  the  2  to 
5-kilometer  altitude  bracket  would  extend  the  usefulness  of  the 
Z  *  200  approximation  to  sotnewhat  higher  rainfall  rates. 

7.0  ALllTODE  DEPENDERCE  OF  THE  EFFECTIVE  REFLECnVITY  OF  SEVERE 

CONVECTIVE  STORMS, 

The  altitude  dependence  of  the  effective  *-eflectlvity  of 
severe  convective  storms  has  been  well  doctssentud  in  the  litera¬ 
ture.  (8,  37,  38,39),  Several  of  the.?*  storms  passed  through  the 
propagation  path  area  during  the  period  of  the  FOPSI  Project  but 
only  one  of  them  has  been  analyzed  in  detail.  This  storm  occurrcsd 
on  28  June  1966  and,  although  not  accompanied  by  extremely  heavy 
surface  rainfall  resulted  in  effective  reflectivities  in  the  ovder 
of  10^  nm  ^m“3.  The  altitude  intervals  and  magnitudes  of  the 
reflectivities  noted  in  connection  with  this  storm  are  shown  in 
Figtires  7,  8,  and  9.  Because  the  common  volume  samples  were  selected 
at  random  rather  than  as  a  result  of  probing  the  storm  for  reflectivity, 
the  actual  maximum  reflectivities  existing  in  the  cells  may  have 
been  somewhat  higher  than  indicated.  During  this  storm  the  fading 
rate  (Figure  10),  the  refractive  index  profiles  (Figures  11,  12,  end 
13)  and  the  photographs  of  the  PPI  scope  of  the  WSR  57  radar  at 
NAFEC  (Figures  14,  15,  and  16)  clearly  indicate  precipitation  scat¬ 
tering  as  the  mode  of  propagation.  From  the  time  the  front  entered 
the  U.S.W.B.  rain  gauge  system  at  1100  hours  until  its  departure 
at  2100  hours,  the  maximum  rainfall  rate  derived  from  any  of  the 
rain  gauge  samples  was  60  raa/hr.  and  from  the  entire  network  of 
seven  gauges  there  were  only  seven  15-mlnute  samples  that  indicated 
surface  rainfall  rates  in  excess  of  20  nm/hr.  The  high  reflectivltes 
which  persisted  for  some  time  were  probably  due  to  hall  aloft, 
although  very  little  of  it  reached  the  ground.  The  "fingered" 
appearance  of  the  storm  cells  on  the  PPI  scope  photographs  and  the 
extremely  high  "radar  tops"  as  noted  by  the  NAFEC  weathei*  radar 
operator  (over  17  kilometers)  indicated  the  presence  of  large  hail¬ 
stones  in  this  storm, ^2,  43,  44,  45) 

Since  the  convective  storms  are  the  principal  contributors 
to  the  high  reflectivity  portion  of  the  probability  distribution 
function  for  precipitation  scattering,  and  exhibit  the  greatest 
disparity  between  reflectivity  aloft  and  nearby  surface  rainfall 
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rates,  they  inu3t  be  taken  into  account  in  the  consideration  of 
possible  interference  to  high-reliability  services.  It  is  possible 
that  the  effect  of  th->se  storms  can  be  accomodated  by  a  suitable 
modification  of  an  eaaplrical  Z-R  relationship  or  by  "blocking"  the 
troublesome  altitude  interval  to  common  volumes  for  allocation 
purposes  or  by  introducing  some  other  altitude-discriminant  meteo¬ 
rological  parameter  such  as  the  "melting  level"  or  QOC.  isotherm. 
The  some\^t  complicated  chart  exhibited  in  Figure  17  represents  an 
attempt  in  this  direction. 

8.0  MmiFICATION  OF  THE  Z-R  RELATIONSHIP  TO  ACCOMODATE  THE  HIGH 
REFLECTIVITIES  OF  C(»IVECTIVE  STORMS . 


As  previously  pointed  out,  the  elimination  of  common  volumes 
in  the  two  to  five  kilometer  altitude  interval  in  the  New  Jersey  area 
would  extend  the  usefulness  of  the  Z  =  200  R^"^  approximation  to 
the  point  where  an  acceptable  level  of  precipitation  scattered 
interference  could  be  predicted  from  the  surface  rainfall  rate 
probability  distributions.  However,  when  the  interference  proba¬ 
bility  prediction  was  .extended  to  p  .0005,  without  regard  to 
coamon  volume  altitude,  the  application  of  the  Z  =  200  R^*^  approxi¬ 
mation  resulted  in  estimates  of  Z  from  R(p)  approximately  10  dB 
below  those  indicated  by  the  New  Jersey  measurements. 

By  adjusting  the  rainfall  rate  distribution  function  to  fit 
the  measured  Z  distribution  function  at  points  of  equal  probability 
it  was  possible  to  arrive  at  a  Z(p)  -  R(p)  relationship  which  more 
nearly  describes  that  of  the  New  Jersey  data.  This  adjustment  was 
accomplished  by  the  least  squares  fit  of  the  regression  line 

log  Z(q)  =  log  a  +  b  log  R(q)  (8-1) 

where  q  «  p(R)  for  1.0  <  R  <  105  naa/hr. 

The  relationship  obtained  in  this  manner  was 

Z(p)  =  127.7  R(p)^'^-^  (8-2) 


and  is  shown  in  Figure  18  along  with  the  regression  line  represented 
by  Z  =  200  R^'”  and  some  recent  Z-R  relationships  from  the  dense 
Illinois  rainfall  suirvey  system. The  trend  of  the  Z  values  in 
the  neighborhood  of  10^  indicates  that  even  this  new  relationship 
would  probably  result  in  the  underestimation  of  Z  (effective)  for 
the  high  moisture  content  zones  of  severe  convective  storms.  However, 
this  approximation,  when  applied  to  the  New  Jersey  data,  results  in 
a  standard  deviation  for  Z(p)  of  less  t^n  1.6  dB  and  appears  to  be 
more  representative  than  is  Z  =  200  R^'®  for  the  actual  relationship 
between  Z(p)  and  R(p)'for  convective  storms. 
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9.0  THE  EFFECT  OF  RAINFALL  INTEGRATION  INTERVAL  UPOS  THE  PROBABILITY 

DISTRIBUTION  F’jUCTIO.i. 

Because  the  relationship  Z(p)  =  127,7  S{p)^*^^  was  derived  by 
considering  15-niinute  rainfall  accumulations  rdther  then  the  hourly 
accumulations  normally  aAmllable  in  most  areas,  there  arises  a  prob* 
lem  in  connection  with  the  application  of  that  relationship  or  any  , 

such  relationship  in  a  practical  allocations  plan.  (47,  48,  49,  59) 

There  appears  to  be  an  additional  requirement  for 

(1)  A  statistical  link  between  the  probability  distribution 
functions  for  various  integration  periods,  or 

(2)  A  showing  that  the  difference  in  probability  distribution 
functions  is  insignificant  within  the  limits  of  the 
extrapolation  of  the  data. 

Basically,  the  determination  of  a  rainfall  rate  from  the  accumulation 
for  any  given  period  amounts  to  obtaining  the  mean  or  average  value 
of  the  Instantaneous  rainfall  in  a  sample  size  corresponding  to  the 
length  of  the  period.  The  distribution  of  rainfall  rates  computed 
in  this  manner  is  actually  the  sampling  distribution  of  the  meane 
(or  averages).  Since  the  population  is  infinite  for  all  practical 
purposes  and  the  number  of  samples  is  large,  the  sampling  distribu¬ 
tion  of  the  means  must  be  approximately  noniial  regardless  of  the 
size  of  the  samples  (length  of  the  accuniulatlon  period)  and  the  dis¬ 
tribution  of  the  population  (instantaneous  rainfall  rates)  itself.* 
since  one  has  almost  no  chance  of  determining  the  instantaneous 
rainfall  rates  either  directly  or  by  the  manipulation  of  rainfall 
accumulations,  and  is  almost  certain  to  encounter  difficulties  in 
determining  the  distribution  function  for  the  high  rainfall  rates 
which  occur  in  very  short  time  Intervals  in  temperate  climates,  one 
is  forced  to  the  alternative  (2),  above. 

To  evaluate  the  effect  of  sampling  interval  differences  upon 
the  rainfall  distribution  function,  a  cumulative  probability  dis¬ 
tribution  of  the  hourly  accumulations  of  the  U.  S.  W.  6.  rain  gauges 
at  Cape  I4ay,  Freehold,  Glassboro,  Hightstown,  Lumberton,  Marlton, 
and  NAFEC  was  generated  for  the  same  time  periods  covered  by  the 
probability  distribution  derived  from  the  15-minute  accumulations. 

An  adjustment  of  this  distribution  function  at  points  of  equal 

probability  by  a  least  squares  fit  of  the  regression  line  (8-1)  ■ 

for  p(R)  1.0 <  R  <  Max.  mm/hr.  resulted  in  the  relationship:  j 

Z(p)  =  99.4  R(p) (9-1) 


*  /  Special  case  of  central  limit  theorem  of  probability  theory. 
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Th«  application  of  the  approximation  developed  from  the 
15-iaiaiite  rain  gauge  data  Z(p)  »  R(p)^*^^  to  the  dis” 

tribut-or  function  of  the  hourly  data  increases  the  stacsiard 
deviatici.  for  Z(p)  fros"  1.57  to  2.38  dB  for  the  Nev  Jersey  data. 

Based  upon  this  i-riterion  alone,  one  might  be  persuaded  to  use  the 
approximation  resulting  frosa  the  least  squares  fit  of  the  regression 
line  of  the  hourly  rain  gau^e  data  In  spite  of  tlie  fact  that  the 
distribution  contains  only  one  quarcer  the  number  of  samples  found 
in  the  distribution  for  the  shorter  time  period  and  the  maximum  rain¬ 
fall  rate  is  reduced  from  over  i.00  icn/hr.  to  approximately  40  nan/hr. 
Manipulation  of  the  time  increments  of  the  rainfall  accunmlations 
of  the  seven  U.  S.  W.  B.  rain  gauges  used  in  the  POPSI  Project  does 
not  represent  an  unbiased  test  of  either  the  relationship  of  (8-2) 
or  that  of  (9-1),  A  more  significant  test  vas  laade  by  considering 
the  probability  distributions  for  the  hourly  precipitation  accumula¬ 
tions  obtained  from  the  Local  Climatological  Data  for  the  U.  S.  W.  B. 
stations  at  Atlantic  City,  Philadelphia,  Newark,  and  Trenton,  during 
the  period  covered  by  the  POPSJ  data.  These  precipitation  data  are 
in  the  form  readily  available  in  most  areas.  In  Figure  19  the  proba¬ 
bility  distributions  of  the  hourly  rain  gauge  accumulations  are  com¬ 
pared  with  those  derived  frosi  Che  15-minute  accumulations.  The  least 
squares  fit  regression  lines  of  log  Z(p)  upon  log  R(p)  are  shown  in 
Figure  20  and  summarized  in  Table  II  which  indicates  the  "goodness 
of  fit"  of  each  relaticoahlp  to  the  data  from  which  it  was  derived. 
Table  III  summarises  the  results  of  testing  the  various  Z(p)  -  R(p) 
approximations  on  rain  gauge  data  other  than  that  from  which  each 
relationship  was  derived.  Although  the  differences  in  the  standard 


errors  of  estimate  of  each  approximation  are  not  really  significant 


the  teat  on  Independent  data  Indicates  that  the  Z(p)  *  127.7  R(p) 


approximation  Is  the  best  within  the  limits  of  the  extrapolation  of 
the  New  Jersey  data. 


10,  CONCLUSION 

The  detailed  analyses  of  the  POPSI  Project  precipitation 
data  have  resulted  in  a  useable  Z-R  relationship  based  upon 
measurements  of  the  scattered  elecjtro^gnetic  field.  By  using 
the  relationship  Z(p)  =  127.7  R(p)^'^”,  in  conjunction  with  the 
probability  distribution  of  high  surface  rainfall  rates,  the 
probability  of  exceeding  a  critical  effective  radar  reflectivity 
in  a  given  common  volume  can  be  calculated. 


The  problems  of  the  altitude  variability  of  the  effective 
reflectivities,  discrete  scattering  volume  definition  and  identi¬ 
fication  of  the  exact  nature  of  the  scatterers  dictate  against  the 
eiTtrapolation  of  low  rainfall  rate  Z-R  relationships  to  the  high 
rates  which  are  pertinent  to  the  calculation  of  possible  interference 
due  to  precipitation  scattering. 
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